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We describe the theory of few Coulomb-correlated electrons in a magnetic quantum dot formed 
in graphene. While the corresponding nonrelativistic (Schrodinger) problem is well understood, 
a naive generalization to graphene's "relativistic" (Dirac-Weyl) spectrum encounters divergencies 
and is ill-defined. We employ Sucher's projection formalism to overcome these problems. Exact 
diagonalization results for the two-electron quantum dot, i.e., the artificial helium atom in graphene, 
are presented. 
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The recent spectacular progress in preparing and use- 
fully employing indivual carbon monolayers of graphene 
U continues to stimulate much interest across differ- 
ent scientific communities, including material science, ap- 
plied physics, chemistry, condensed matter physics and 
mathematics. Ballistic electronic motion with quantum 
coherence extending over micrometer distances has been 
achieved in several experiments, see, e.g., Ref. Q. The 
low-energy physics close to a single K point can then be 
described by a two-component Dirac-Weyl Hamiltonian 
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H a = v-p <t ■ [p — -A 
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suggesting an easily accessible condensed-matter realiza- 
tion of relativistic quantum mechanics. In Eq. {T]), er 
denotes the vector of the first two Pauli matrices for the 
"isospin" encoding the two sublattices, the Fermi velocity 
is vp m 10 6 m/sec, and we include a static vector poten- 
tial A{r) describing (possibly inhomogeneous) magnetic 
fields. Since graphene's effective fine structure constant 
is a w 1, present interest has also turned to Coulomb 
interaction effects @]. According to recent Monte Carlo 
simulations [j| and analytical arguments [||, sufficiently 
strong interactions may even open a sizeable bulk gap in 
the Dirac fermion spectrum. 

Here we study the properties of few Coulomb- 
correlated electrons confined to a finite-size quantum 
dot formed in graphene. Using electrostatically formed 
quantum dots in semiconductor devices, such "artificial 
atoms" have been intensely studied over the past two 
decades, both experimentally Q and theoretically Q. In 
graphene dots formed by electrostatic gating, however, 
carriers can escape due to the (recently observed @) 
Klein tunneling phenomenon, and at best quasi-bound 
states may appear [Io| . An alternative is to employ 
lithographically defined quantum dots fill ], where de- 
tailed information on ground- and excited-state proper- 
ties has been obtained from transport spectroscopy. Un- 
fortunately, the boundary of lithographically fabricated 
gra phene dots is rather disordered and difficult to model 
fl2| . On the other hand, suitable and realizable inho- 
mogeneous magnetic fields can confine Dirac fermions 



promising to yield tunable and well-defined mag- 
netic graphene dots. 

A more challenging difficulty to theory arises when 
trying to generalize Eq. (|TJ) to a first-quantized many- 
particle description. The first-quantized approach has 
turned out to be very efficient and convenient for the 
case of Schrodinger electrons in semiconductor-based ar- 
tificial atoms 0]. For the "relativistic" graphene case, 
the problem arises from the unboundedness of Eq. ([1]), in 
contrast to the corresponding Schrodinger operator (p — 
- c A) 2 jlm* . While Eq. |T]) can still be used within effec- 
tive single-particle approximations such as the Hartree- 
Fock approach 151, variational schemes [Hj], or density- 
functional theory [lj], the full iV-particlc problem (for 
small N > 1) with Eq. fTJ) for the kinetic part suffers 
from the so-called "Brown- Ravenhall disease" [ill, [Tj|. 
Roughly speaking, the unbounded spectrum allows par- 
ticles to lose arbitrary amounts of energy by transferring 
their energy in (real) scattering events to other particles. 
The resulting divergent density of states prohibits, for ex- 
ample, the direct use of exact diagonalization methods. 
This difficulty of the Dirac equation has been known for 
half a century . To "cure" this "disease" , we follow a 
proposal by Sucher [l9[ and confine the Hilbert space to 
positive energy eigenstates through suitably defined pro- 
jectors, cf. Eq. © below. While we formulate this ap- 
proach for the magnetic dot only, the general concepts re- 
main applicable for almost arbitrary graphene dots. The 
projection method then allows, for instance, to apply nu- 
merical techniques to the relativistic ./V-particle problem. 
In this work, we present exact diagonalization results for 
the many-body energy spectrum of the artificial helium 
atom (N = 2) in graphene. 

Let us first specify the model discussed here (we set h = 
1). In cylinder coordinates, we consider the spherical and 
parabolic magnetic field profile oriented perpendicular to 
the graphene plane (with A in symmetric gauge), 



c x , -smip 

B(r,<p) = -ulr 2 e z , A(r,<p) = -w|r 3 cos^ 
e 4 \ 

'(2) 

The inverse length scale -Jlub tunes the field inhomo- 



geneity. The dimcnsionlcss radial coordinate is g = 
rJws, and energy (e) is measured in units of vf^/lob- 
[Physical units are recovered from \J B{g) /Tcs\a = 
(vy \Jojb 1 26meV) g] Such magnetic profiles can be gener- 
ated with reasonable accuracy using suitable lithograph- 
ically defined ferromagnetic films deposited on top of the 
graphene layer after formation of a protective oxide layer 

[2(ij |. Upper and lower components of eigenspinors [i/v ) 
to Eq. (JTJ) must then differ by one orbital angular momen- 
tum quantum number (to) due to conserved total angular 
momentum [HI ]. With real functions 4> m (Q) and Xm+i(o), 
the radial part of [# ) is oc (<j) m (g),i sgn(e) Xm+i(o)) T , 
where Eq. |T]) yields the radial equations 
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which cannot be solved analytically. We here carry out 
exact diagonalizations, later on including the Coulomb 
interaction, and thus solve Eq. ([3]) numerically. It is con- 
venient to employ the Darwin-Fock states (with integer 

n > and the Laguerre polynomials Ln™') [H 
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as complete orthonormal function set in two dimensions 
(2D) to expand <f) m and Xm+i in Eq. ((3j) . In Eq. (f4|), 
we have included an additional tunable width parameter 
A, which can be optimized [2l[ to reduce the number of 
required basis functions < n < n max when approximat- 
ing the to the desired accuracy. Figure Q] displays 
the resulting eigencnergies as a function of the orbital 
angular momentum to. As expected, the spectrum is 
electron-hole symmetric, and for m > 0, a zero-energy 
level develops. This zero-energy level is non-dispersing 
(precisely like a quantum Hall level), despite of the in- 
homogeneous magnetic field which implies the nontriv- 
ial TO-depcndence of all other energy levels. Note that 
for the corresponding Schrodinger case with a parabolic 
magnetic field, the zero-energy level is absent. 

Next wc consider N interacting electrons in such a 
graphene dot. A naive approach is to consider the first- 
quantized Hamiltonian 
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H = v F J2<Tj-\Pj-Mrj)]+J2 



CtVF 



(5) 



i<3 



where the fine structure constant is a = e 2 /(noVF)- 
For typical substrate materials, the dielectric constant is 
kq w 1.4 to 4.7, resulting in a w 0.6 to 2. We mention in 
passing that the "Wigner molecule" regime [22j seems out 
of reach in graphene dots, since both the kinetic and the 
potential e nerg y show identical scaling when changing 
the density [231 ] - Moreover, we neglect the Zeeman term 
which is very small in graphene [14| . Up to the spin and 




FIG. 1: (Color online) Exact diagonalization results for 
the single-particle eigenenergies e vs orbital angular mo- 
mentum m in a parabolic magnetic quantum dot in 
graphene [see Eq. (j3}]. The zero-energy levels are indi- 
cated as red filled circles, while the other levels are shown 
as black filled diamonds. 



i4T-point indices 2J|, many-body spinors then have 2 
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components. For the related Schrodinger dot (Hs), con- 
finement of electrons is usually achieved by a parabolic 
electrostatic potential Q , and the many-particle descrip- 
tion analogous to Eq. ((5|) simplifies considerably owing 
to the generalized Kohn theorem [25| . According to this 
theorem, Hs = H cm + H rc i separates into two commuting 
parts describing center-of-mass (H cm ) and relative (H Te \) 
motion. Then H cm is just a 2D harmonic oscillator, while 
H IC \ contains all Coulomb interaction effects. In addition, 
-ff r oi conserves angular momentum, as does H cm . Taking 
N = 2 as example, in effect only a ID quantum problem 
for the radial motion of H Ic i remains to be solved. In 
contrast, Eq. |5]) does not benefit from Kohn's theorem, 
and only the total angular momentum remains conserved 
as dictated by isotropy. Therefore, while the additional 
spinor structure already increases the rank of the Hamil- 
tonian matrix in the Dirac-Weyl case by a factor 2 N , the 
rank grows even more severely because neither H cm nor 
a conserved angular momentum of H TC \ can be separated 
off the problem. For N = 2 (graphene helium), wc needed 
to include states up to n max ps 14 to reach sufficient accu- 
racy. In addition, contrary to the Schrodinger problem, 
particles may now exchange relative angular momentum 
Am through the interaction. Owing to the exponential 
decay of Coulomb matrix elements with |Am|, it is suf- 
ficient to take | Am | < 3, yielding an additional factor 
7 N ^ X to the matrix size (for N particles). For N — 2, we 
then need to include 14 2 x 7 x 2 2 = 5488 product basis 
states in total. 

Let us then address the more fundamental difficulty 
arising already for N = 2 when naively using Eq. ([5]). A 
closely related pro blem has been pointed out by Brown 
and Ravcnhall |18| in a relativistic treatment of the he- 
lium atom: the Dirac equation analogous to Eq. (T5|) does 
not possess normalizablc antisymmetric eigenstatcs in 



the two-particle Hilbcrt space. This failure has its origin 
in the unbounded spectrum of the Dirac Hamiltonian, 
which allows for unlimited energy exchange among the 
particles. As a result, the density of two-particle states 
increases with Hilbert space dimension and ultimately 
diverges. This causes, e.g., divergent contributions in 
second-order perturbation theory. In consequence, nei- 
ther the true Dirac equation nor the two-component vari- 
ant {TJ) for graphene allow for naive many-particle gen- 
eralizations such as Eq. ((SJ). To overcome this deficiency, 
Sucher proposed to restrict the (anti-symmetrized) 
product Hilbert space to the positive energy eigenspace 
for each particle by means of a suitable projector A+. 
We thus consider a situation with chemical potential 
fx = + just above zero, where all non-positive energy 
(hole) states up to (and including) zero energy are filled. 
As long as a finite energy gap separates the relevant posi- 
tive energies (which can then be occupied by the N elec- 
trons under consideration) from the filled non-positive 
energy levels, such a projection seems physically sensi- 
ble, at least for interaction strengths not exceeding this 
gap. In natural helium, the mass gap ensures the validity 
of such an approach. On the other hand, in a homoge- 
neous 2D graphene sheet, weakly interacting fermions are 
gapless and Sucher's approach does not apply. For our 
finite-size quantum dot, however, there is a finite gap 
between e = and the lowest positive-energy level, see 
Fig. [T] We therefore consider two additional electrons 
confined in the quantum dot, residing above a filled Dirac 
sea with [i = + . The projector proposed by Sucher is 
expressed as 

A + = A«®Af, = $>£'' 0) ><^'' 0) l, (6) 

where the sum is restricted to strictly positive single- 
particle energies 

E v£j_ > ( for particle j = 1, 2) indexed 
by v € /, sec Fig. [TJ with corresponding eigenspinors 
\ipu )■ With H in Eq. (J5j) , the projected Hamiltonian 
i?D = A + HA + is well-behaved and exhibits a finite den- 
sity of two-particle states which does not increase with 
Hilbert space size. Throughout the experimentally rele- 
vant regime, a < 2, interactions are not strong enough 
to induce a breakdown of this projection approach. 

We have carried out exact diagonalizations of Hd for 
N = 2 using a two-particle product basis of Darwin- 
Fock states (UJ). This basis has the advantage of allow- 
ing to analytically express the matrix elements of the 
two-particle interaction operator in Eq. ([5]) in terms of 
finite sums and products, i.e., no quadratures nor trun- 
cations of infinite sums are necessary; for the correspond- 
ing (lengthy) expressions, see Ref. [26|. Their numerical 
evaluation involves taking small differences of huge num- 
bers, the latter increasing as n max ! with the number of 
kept Darwin-Fock states. We employed algorithms al- 
lowing for number manipulations of arbitrary precision 
and used 30-digits accuracy. The resulting energy spec- 
trum of artificial helium in a magnetic graphene quantum 
dot is shown in Fig. [2] for conserved total orbital angular 











: , '-t^ 1 ;:- 1 ^- 






M=-2 






M=-l 






M=0 






M=l 




<r i , i 


1 



0.5 1 1.5 

a 

FIG. 2: (Color online) Exact diagonalization results for 
the energy spectrum of graphene artificial helium vs fine 
structure constant a for different total angular momenta 
M. States with M = — 2 are shown as black solid, M = 
— 1 as blue dashed, M = as green dot-dashed, and 
M = 1 as red dotted curves. 

momenta M = mi + m 2 = —2,-1,0,1. These values 
include the ground state for a < 2. All levels rise when 
increasing the (repulsive) interaction strength. This also 
holds true for the hole states (not displayed in Fig. [2]), 
which, however, cross e = only for a>2. For a < 2, 
interaction matrix elements indeed remain much smaller 
than the energy difference Ae ~ 2.98165 between the 
lowest two-particle state (for M = — 2 at a = 0) and 
the zero-energy level, a posteriori justifying Sucher's ap- 
proach here. Figure [5] reveals that states with larger total 
angular momentum M, or higher excited states, tend to 
increase less in energy with a as compared to the M = — 2 
ground state. This is a consequence of the larger spa- 
tial extent of excited-state wavefunctions, with a reduced 
Coulomb repulsion between the electrons. Particularly 
striking is the shallow increase of the lowest M = 1 en- 
ergy level, which even becomes lower in energy than the 
M = —2 level for a > 1.6. Approximating this level at 
a = by Darwin-Fock levels, one of the two particles 
is seen to have m — 2 for the lower spinor component, 
cf. Eq. ©, causing a significantly larger spread of this 
part of the wavefunction compared to M = —2,-1,0. 
Figure [2] also reveals nontrivial spin physics. In the pres- 
ence of interactions (a > 0), doubly degenerate noninter- 
acting (a = 0) energy levels will split into a spin-triplet 
(S = 1) state of lower energy and a spin-singlet (S = 0) 
state of higher energy, in accordance with Hund's rule. 
The triplet states are approximately (see below) Zeeman- 
degenerate. Singly degenerate a = levels, such as the 
M = —2 ground state (for small a), are 5 = states and 
remain unsplit for a > 0. Thus we expect singlet-triplet 
ground-state spin transitions to occur within < a < 2, 
as the one seen in Fig. HataK 1.6. 

Finally, we remark on optical transitions between the 
many-body energy levels in Fig. [5] For the electro- 
statically defined parabolic Schrodingcr quantum dot, 
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the generalized Kohn theorem implies that Coulomb in- 
teractions can never affect optical transitions because 
the dipole operator J2j=i r j &c ^ s exclusively on the 
eigenspace of H cm . Therefore optical spectra just reflect 
the harmonic excitations of the center-of-mass motion 
Q. However, in our magnetic graphene dot, Kohn's the- 
orem is ineffective and optical transitions between dif- 
ferent many-body levels in Fig. [5] are possible, thereby 
allowing to optically probe interaction physics. Note 
that magnetic fields are usually assumed homogeneous 
such that photons cannot change the total spin S of the 
charged many-particle system in electrical dipole transi- 
tions. While this would prohibit all transitions between 
states with different S, the inhomogeneous magnetic field 
here (slightly) mixes the S z = components of S = 
and S = 1 levels. We estimate the amount of this mixing 
by the variation of the Zeeman energy across the spatial 
extent of the wavefunction compared to the level sepa- 
rations of Hy>. As a first estimate, compare the Zeeman 
energy Az at the maximum of the charge-density distri- 
bution with the orbital (Landau) energy at this point, 
A Z /A L = g/j, B B/V2ecB ~ 10" V B / Tcsla - While this 
is small, the Zeeman energy variations can easily exceed 



orbital level separations near spin singlet-triplet degen- 
eracies, e.g., for a <C 1 or close to level crossings in Fig. [2j 
resulting in a strong spin mixing. The corresponding 
transitions are then optically allowed. 

To conclude, we have presented the theory of few in- 
teracting electrons in a (magnetically confined) graphene 
quantum dot. The low-energy Dirac-Weyl spectrum of 
graphene suggests that one can realize relativistic artifi- 
cial atoms in this setting. While a naive formulation en- 
counters conceptual difficulties related to the unbounded- 
ness of the Dirac-Weyl Hamiltonian, by virtue of Sucher's 
projection operator approach, a consistent and accurate 
theory can be given. We have presented exact diago- 
nalization results for the energy spectra of artificial he- 
lium, where we predict singlet-triplet ground-state spin 
transitions to occur for a < 2. Moreover, the reported 
many-body levels can be experimentally probed by opti- 
cal spectroscopy. — We thank H. Siedcntop for drawing 
our attention to Refs. and and acknowledge dis- 
cussions with A. Dc Martino, K. Ensslin, and T. Hcinzel. 
This work was supported by the DFG (SFB TR/12) and 
by the ESF network INSTANS. 
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